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We use numerical simulations to show how noninteracting hard particles binding to a deformable 
elastic shell may self-assemble into a variety of linear patterns. This is a result of the nontrivial 
elastic response to deformations of shells. The morphology of the patterns can be controlled by the 
mechanical properties of the surface, and can be fine-tuned by varying the binding energy of the 
particles. We also repeat our calculations for a fully flexible chain and find that the chain confor- 
mations follow patterns similar to those formed by the nanoparticles under analogous conditions. 
We propose a simple way of understanding and sorting the different structures and relate it to the 
underlying shape transition of the shell. Finally, we discuss the implications of our results. 



I. INTRODUCTION 

Spontaneous assembly of components into large or- 
dered aggregates is a ubiquitous phenomenon in nature, 
and is observed across all length scales. Aggregation of 
proteins into functional nanomachines pQ, formation of 
viral capsids [21 13] , packing of phospholipids into biologi- 
cal membranes 1) and assembly of colloidal particles into 
macroscopic photonic crystals [5] are just a few examples 
of such fundamental processes. The human body itself 
develops as a hierarchical self-assembly of units of differ- 
ent sizes. Understanding the physical mechanisms driv- 
ing self-assembly of biological or artificial components 
will shed light on critical biological processes and simul- 
taneously holds promise for the development of materials 
with novel functional, mechanical, and optical properties. 

The use of fluid interfaces has revolutionized the field 
of self-assembly [SI E] by providing a universal mech- 
anism through which arbitrary building blocks can be 
driven close to each other. As a result the process of self- 
assembly, once specific to molecular systems, has been ex- 
tended to components of much larger size [B] . The physics 
behind this mechanism is well known [8l [9] , and its man- 
ifestations are a matter of common experience, as seen 
for instance in the clustering of cereal on milk surface. 
Local deformations in the interface profile, caused by the 
presence of floating objects, induce long-range capillary 
forces. These are the result of the minimization of the 
total surface area of the interface that is regulated by its 
surface tension. By manipulating the surface tension of 
the interface, and by tailoring the chemistry of the build- 
ing blocks, millimeter-size objects, microchips and mi- 
crocomponents have been successfully self-assembled [TDJ- 

In this work we try to understand how elastic sur- 
faces, which have quite different mechanical properties 
from their fluid counterparts, can be used to drive self- 
assembly of nanoparticles. Unlike fluid interfaces, whose 
molecules can freely diffuse across their area, the ele- 
ments of elastic surfaces are tethered together, and re- 
spond to small deformation in a spring-like fashion. The 
resulting macroscopic behavior, regulated by bending 



and stretching energies, is characterized by strong non- 
linearities |13| . The mechanical properties of macro- 
scopic elastic sheets have recently been the subject of 
intense scrutiny 11 19 . However, there are several mi- 
croscopic artificial and naturally occurring examples of 
such materials, including graphite-oxide sheets [201 I2T] . 
graphene sheets [221 [23], cross polymerized biological 
membranes [23], buckypaper [25rl2T] . the spectrin-actin 
network forming the cytoskeleton of Red Blood Cells [2H1 
125] . and very recently they have been fabricated using 
close-packed nanoparticle arrays [50] . 

It is easy to show [T51 HO] that the ratio between the 
stretching and bending energies for an arbitrary deforma- 
tion of amplitude h on a flat elastic sheet of thickness t 
scales as E s /Eb ~ (h/t) 2 . Therefore, for sufficiently thin 
sheets, bending is the preferred mode of deformation. 
This has a profound effect on the way elastic surfaces 
respond to deformations as the only stretch-free defor- 
mation involves uniaxial bending. Skin wrinkling under 
applied stress [El [17] and stress focusing via d-cone for- 
mation of crumpled paper [19] are two beautiful examples 
of this phenomenon. We have recently shown how the re- 
sponse to deformations of elastic nanotubes can alter the 
elastic properties of a flexible filament binding to it [31j . 
and drive self assembly of colloidal particles into helical 
and axial lines [52] , 

In this paper we show how the elastic response to defor- 
mations of a spherical shell can be used to self-assemble 
colloidal particles in a variety of patterns that only de- 
pend on the mechanical properties of the system and the 
amount of deformation of the surface. We also show how 
a fully flexible polymer bound to the shell will sponta- 
neously arrange to conform to similar patterns observed 
for the nanoparticles. What makes this problem very in- 
teresting is that unlike planar or cylindrical geometries, 
for which there is a clear solution to the stretch-free de- 
formation problem, any deformation of a spherical shell 
will necessarily involve stretching of the surface. The 
shape of the deformation minimizing the stretching en- 
ergy in this case is therefore not immediately obvious. 
Furthermore, a buckling transition from the spherical to 
a faceted icosahedral shape is known to take place at 
large stretching energies |33| . and finally for sufficiently 
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small shells (or large deformations) the ratio between 
bending and stretching energy becomes independent of 
h, and only depends on the radius of the sphere R and 
its thickness t, E s /E b ~ (R/t) 2 [13]. 



II. MODEL 

We model the elastic surface via a standard fishnet 
network representation [31] . Each node of the network 
is placed to conform to the symmetry of an icosadelta- 
hedron. In such surfaces, like viral shells, all but 12 
nodes have a regular triangulation with six neighbors; 
12 five-fold disclinations are also present as required 
by Euler's theorem relating number of faces, edges and 
faces on a spherical triangulation. The number of nodes 
Nk describing the surface is then related to the loca- 
tion of the five disclinations, and satisfies the constraint 
Nk = 10(n 2 + nm + to 2 ) + 2 [351. Here n and to indicate 
that one must move n nodes along the row of neighboring 
bonds on the sphere, and then after a turn of 120°, move 
for m extra steps. 

We studied two different sphere sizes; the smaller one 
has Nk — 6752 nodes and symmetry described by n = 15 
and m = 15, the larger one contains Nk — 10832 nodes 
with n = 19 and to = 19. 

To impose surface self-avoidance we place hard beads 
at each node of the mesh. Any two surface beads interact 
via a purely repulsive truncated and shifted Lennard- 
Jones potential 
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where r is the distance between the centers of two beads, 
c is their diameter, and e = 100fc B T. 

We enforce the surface fixed connectivity by linking 
every bead on the surface to its first neighbors via a har- 
monic spring potential 



Ustretching = K s (r — V b) 



(2) 



Here K s is the spring constant and r is the distance be- 
tween two neighboring beads, tb — 1.23c is the equilib- 
rium bond length, and it is sufficiently short to prevent 
overlap between any two triangles on the surface even for 
moderate values of K s . 

The bending rigidity of the elastic surface is modeled 
by a dihedral potential between adjacent triangles on the 
mesh: 
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where <f> is the dihedral angle between opposite vertices of 
any two triangles sharing an edge and Kb is the bending 
constant. 

Colloidal particles of diameter c c = 10c are described 
via the repulsive truncated-shifted Lennard-Jones po- 
tential introduced in Eq. [I] with c — > c c . Finally the 



polymer is constructed as a "pearl necklace" with N £ 
[20, 45] monomers of diameter also c m = 10c. Neigh- 
boring monomers are connected by harmonic springs 
as in Eq. [2] with the equilibrium bond length rj; = 
1.18c m and spring constant of 120fceT/c 2 . Polymer self- 
avoidance is again enforced via the repulsive truncated- 
shifted Lennard-jones potential introduced in Eq. [TJwith 
c -> c m . 

The generic binding of the polymer monomers (and the 
colloids) to the surface is described by a Morse potential: 
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where r is the center-to-center distance between a 
monomer and a surface-bead, Tmb is bead-monomer con- 
tact distance tmb = 5.5c and Dq is the binding energy. 
The interaction cutoff is set at 10c and 7 = 1.25/c. 

We used the LAMMPS molecular dynamics pack- 
age [36] with a Nose/Hoover thermostat in the NVT 
ensemble to study the statistical behavior of the system. 
The timestep size was set to dt = 0.002to (to is the di- 
mensionless time) and each simulation was run for a min- 
imum of 5 TO 6 steps. The radii of the undeformed spheri- 
cal shells are R = 29.05c for N k = 6752, and R = 34.16c 
for N k = 10832. 



III. RESULTS 

To understand how the configurational properties of 
the binding polymer and colloids are related to the elas- 
tic properties of the templating surface, we performed a 
series of simulations for many values of K s , Dq and Kb, 
and for different polymer lengths and number of colloidal 
particles. We find that a convenient way of representing 
our data is via the dimensionless parameter known as 
the Foppl-Von Karman number, defined as 7 = YR 2 /k, 
where Y is the Young modulus of the shell and k is the 
bending rigidity as defined in the continuum theory of 
elasticity [T3J. We begin our analysis by studying the 
buckling transition of the shell as a function of 7 in the 
absence of binding agents. This will give us critical infor- 
mation about how to relate the shape of the templating 
surface (the elastic shell) and its elastic properties. To 
identify the buckling transition and match the numerical 
parameters of our model with 7 we follow the analysis 
carried out in [53] . Fig. [T[\ shows the results of our sim- 
ulations. The surface asphericity A, defined as 
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is plotted againts 7 = (4/3)K s R 2 /K b . The buckling 
transition, for which the spherical shape begins to facet, 
begins for values of 7 > 10 2 . This result is in good agree- 
ment with that computed in |33j , and represents a good 



3 



test of our numerical model. We next add colloidal parti- 
cles to the system and observe their self-assembly on the 
surface of the spherical shell. Depending on the value of 
7, particles arrange in patterns that minimize the elastic 
energy of the shell. Figure |Tj3 shows the resulting pat- 
terns as a function of 7 for different values of indentation 
h, and a constant number of colloidal particles N — 35. A 
convenient way of extracting the indentation is obtained 
by computing (A p — A) 1 / 2 , where A p is the asphericity 
of the shell in the presence of the particles and A is the 
asphericity without them. 

The general feature of the diagram is that self-assembly 
occurs for only a relatively narrow range of particle defor- 
mations. When (A p — A) 1 / 2 is too small, i.e. the surface 
is basically unaffected by the presence of the particles, we 
find no aggregation; this is the result one should indeed 
expect when placing N noninteracting hard particles on 
the surface of a rigid sphere (al low densities). We in- 
dicate this phase as the gas phase. When [A p — A) 1 / 2 
is too large and significant deformations are induced by 
the binding particles, we find that the system becomes 
kinetically trapped (at least within the timescales con- 
sidered in this study) into metastable states. In fact, re- 
peating the simulations under the same conditions leads 
consistently to different and not well defined aggrega- 
tion patterns. We indicate this phase as the arrested 
phase. For intermediate deformations, particles consis- 
tently self-assemble in a variety of patterns whose fea- 
tures are clearly related to the mechanical properties of 
the shell via 7. 

For small values of 7 we find that particles aggregate 
isotropically to form a two dimensional crystal on top of 
the sphere. The presence of these 2d crystals tends to 
flatten the surface underneath it. As a consequence, if 
a sufficiently large number of particles is added to the 
system we find that the side length of the crystal is lim- 
ited by the shell diameter, and extra particles begin to 
aggregate on its opposite side (Fig. [TJ3-I) . As 7 becomes 
larger, at low surface coverage, particles become localized 
over the twelve disclinations, and as N increases linear 
aggregates initially grow by linking the five fold disclina- 
tions on the sphere and finally form a linear network with 
3-line joints winding around the disclinations in the shell 
(Fig. [TJ3-II) . As the templating surface begins to facet, 
each segment of the network straightens revealing a clear 
patterns following the seams that a pentagonal tiling of 
the sphere would generate (a dodecahedron, Fig. [Tf3- 
III). It should be noticed that the previous two phases 
are topologically equivalent. The only difference is the 
presence of extra particles sitting on the twelve disclina- 
tions in region II of the diagram. Increasing the surface 
coverage in region III results in thickening of the width 
of the dodecahedral pattern by formation of parallel and 
adjacent secondary lines of particles. 

At even larger values of 7, once the shell is well faceted, 
particles arrange into a linear and non-connected aggre- 
gate that smoothly winds around and away from the 
twelve disclinations (Fig. [TJ3-IV). This pattern is remi- 



niscent of that of the baseball or tennis ball seam, the dif- 
ference being that in our case the geometry of the aggre- 
gate is dictated by the presence of 12 topological defects, 
while in the baseball the seam winds around the loca- 
tion of the four s = 1/2 disclinations that a thin nematic 
liquid crystal texture would generate on a sphere [37|,l3"8]. 

Finally, for the largest values of 7 we find that particles 
arrange into straighter but shorter aggregates (rods) that 
localize into ten distinct regions of the shell. The length 
of the rods grows with the surface coverage until a critical 
size which depends on the size of the sphere and equals 
roughly the distance between two disclinations located at 
the opposite vertices of two of the icosahedral triangles 
that share one edge (Fig. [l^-V). Further increase of 
iV" results in the formation of multiple rods per region 
which align parallel to each other. Figure [2] shows the 
explicit dependence of the pattern as a function of surface 
coverage in region II of the phase diagram. 

To gain insight into the physical origin of the different 
patterns formed by the particles, we measured the strain 
and bending energy at each node of the shell for differ- 
ent values of 7 in the absence of the particles. Figure [3] 
shows the energy map of the two contributions for small, 
intermediate and large values of 7. The results are quite 
revealing and provide a simple framework from which the 
patterns can be understood. Particles aggregates align 
to follow the low bending and stretching energy regions 
on the shell. The formation of isotropic aggregates and 
the presence of the particles on top of the disclinations 
for small values of 7 suggests that bending energy plays 
an important role in determining the pattern in region 1 
and 2 of the diagram. Regions III, IV and V are instead 
completely dominated by the stretching energy which is 
driving the shape transition of the shell. 

Of particular interest are the two phases that occur 
for very large values of 7 (region IV and V). This con- 
dition that can be obtained either by increasing the ra- 
dius of the sphere R or by significantly altering the rel- 
ative weight of stretching and bending energies in favor 
of the former. In this regime particles can easily bend 
the surface in regions that are far away from the icosa- 
hedral vertices, yet the only bending deformations that 
will not induce stretching energy are those that involve 
bending around a single radius of curvature. We believe 
that the transition from the dodecahedral arrangement 
to the smooth closed loops in Fig. [TJ3 region IV is due to 
the large stretching energy cost that would be associated 
with the formation of either sharp corners or three-lines 
joints. This constraint become so severe in region V of 
the same figure, that the linear aggregates break into 
shorter pieces of roughly equal length that straighten by 
flattening the edges shared by any two triangles of the 
underlying icosahedral geometry. 

Interestingly, we find that almost all phases depicted 
in Fig. [l] can also be acquired by a fully flexible chain 
binding to the shell for analogous degrees of deforma- 
tions and comparable number of monomers. The only 
differences are due to the connectivity constraints on the 
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chain. Region II and region V are obviously impossible 
to achieve with a chain; nevertheless, the difference in re- 
gion V is minimal as short rods are in this case replaced 
by a poly-line with segments having the same length of 
the rods formed by the colloids. In region II the polymer 
traces a simple path connecting the disclinations. Fig- 
ure [4] shows snapshots of the different phases for the fully 
flexible chain. Given these similarities between the be- 
havior of the polymer and that of the colloidal particles, 
it is clear that the mechanism regulating the conforma- 
tional changes of the chain on a flexible shell is identical 
to that driving the self-assembly of the colloidal particles 
on the same surface. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper we detail how deformable elastic surfaces 
can be used to mediate self-assembly of otherwise non- 
interacting colloidal particles, and/or alter the confor- 
mational properties of a fully flexible chain bound to it. 
We find that the structure of the aggregates (the con- 
formation of the polymer) can be understood in terms 
of the mechanical properties of the templating shell via 
the Foppl-Von Karman number 7. Specifically, we have 
shown how there are two distinct regimes: one domi- 
nated by the shell bending energy, and the other by its 
stretching energy. In the former case the shell acquires an 
overall spherical shape and particles localize on top of the 
disclinations and organize to link small bending energy 
regions. In the latter case the shell facets into an icosa- 
hedron with vertices located where the twelve five-fold 
disclinations reside; here particles follow low stretching 



pathways across the shell. 

Crucially, the underlying shell's shape transition de- 
termines the role of the disclinations in the self-assembly 
process. The twelve five folded points attract the binding 
colloidal particles for small values of 7 and repel them in 
the other regime. This result can be of great importance 
for controlling the functionalization of mesoparticles such 
as for instance colloidosomes [35] ■ In fact, we have shown 
how a small number of particles can be localized around 
the twelve disclinations in the bending dominated regime, 
and around the vertices of a dodecahedron (dual to the 
icosahedron) in the stretching dominated regime. 

Our results can be considered as another example or 
an extension of the ideas discussed by Nelson [37l 00] re- 
garding the functionalizability of the four s = 1/2 discli- 
nations in a nematic liquid crystal texture covering a col- 
loidal particle. One of the main differences, apart from 
the overall symmetry of the problem, is that in our case 
the shape of the template is allowed to change in response 
to the elastic strains induced by the presence of the de- 
fects, and that self-assembly (or conformational changes 
in a flexible chain) is driven by the elastic response to 
deformation of the shell. It would be interesting to study 
how particles self-assemble over a deformable shell under 
the overall tetrahedral symmetry provided by a nematic 
liquid crystal texture. 
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FIG. 1: (a) Plot of the asphericity A versus Foppl-Von Karman number 7 indicating the buckling transition in our model of a 
spherical elastic shell in the absence of colloidal particles, (b) Phase diagram indicating how the different aggregates formed 
by the colloidal particles depend on the mechanical properties of the shell (7) and the degree of indentation measured in terms 
of the particle-induced asphericity (A p — A) 1 ^ 2 /N computed using N = 35 particles. The five different phases are indicated 
with Roman numerals, (c) Snapshots from our simulations of the phases indicated in the phase diagram. From top to bottom 
7 = 5.6, 7 = 37.5, 7 = 75, 7 = 225 and 7 = 4500. For the sake of clarity, the particles defining the shell are depicted with a 
larger volume and the regions arund the disclinations are depicted with a lighter color. 



7 





FIG. 3: Bending energy map (TOP) and stretching energy map (BOTTOM) of a spherical shell as a function of 7. Different 
shades indicate the relative strengths as indicated in the color map on the side. 




FIG. 4: Snapshots showing the different conformations adopted by a fully flexible chain when binding to a deformable elastic 
shell. From left to right 7 = 5.6, 7 = 37.5, 7 = f 12.5, 7 = 140 and 7 = 2700. 



